Global microscopic calculations of ground-state spin and parity for odd-mass nuclei 
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Systematic calculations of ground-state spin and parity of odd-mass nuclei have been performed 
within the Hartree-Fock-BCS (HFBCS) approach and the Finite- Range Droplet Model for nuclei for 
which experimental data are available. The unpaired nucleon has been treated perturbatively, and 
axial and left-right reflection symmetries have been assumed. As for the HFBCS approach, three 
different Skyrme forces have been used in the particle-hole channel, whereas the particle-particle 
matrix elements have been approximated by a seniority force. The calculations have been done for 
the 621 nuclei for which the Nubase 2003 data set give assignments of spin and parity with strong 
arguments. The agreement of both spin and parity in the self-consistent model reaches about 80% 
for spherical nuclei, and about 40% for well-deformed nuclei regardless of the Skyrme force used. 
As for the macroscopic-microscopic approach, the agreement for spherical nuclei is about 90% and 
about 40% for well-deformed nuclei, with different sets of spherical and deformed nuclei found in 
each model. 

PACS numbers: 21.10.Hw,21.60.Jz 



To describe the rich variety of elementary modes of nu- 
clear excitations and build up some confidence in the ap- 
proaches used for that purpose, it is important to assess in a 
quantitative systematic way the quality of the reproduction 
of experimental static properties. We choose here to com- 
pare the results of Skyrmc-Hartree-Fock-BCS (hereafter 
referred to as HFBCS) calculations of ground-state (GS) 
spin and parity with experimental values assigned with 
strong arguments available in the Nubase 2003 database [l[ 
for 621 odd-Z and odd-N nuclei. We undertake such a global 
study with three different effective interactions, namely the 
SIII 0, SkM* [J and SLy4 g] parameterization of the 
Skyrme effective nucleon-nuclcon interaction in the mean- 
field channel, together with a seniority force in the pairing 
channel. The results arc also compared with those obtained 
previously in the Finite-Range Droplet Model (FRDM) Q 

We would like to emphasize the global character of the 
present work. Even though some local studies (in limited 
mass regions) of the GS spectroscopic properties of odd- 
mass nuclei using such phcnomcnological interactions have 
been carried out over the years (see, e.g., Refs. 0,0), we 
want to check as completely as possible the relevance of the 
obtained results by limiting ourselves to two observables, 
the spin and the parity of the ground state. The present 
paper is a brief report on these first results. More details 
on the approximations and technical methods will be given 
in a forthcoming study of exotic neutron-rich nuclei. 

Our self-consistent mean-field approach including pair- 
ing correlations relies on the Hartree-Fock and BCS ap- 
proximations. They have been implemented in earlier fis- 
sion studies for even-even heavy and light nuclei [l(| 
with the Skyrme interaction in its SkM* parameterization 
which successfully described potential-energy surfaces in 
terms of constraints on deformation. As for the pairing 
residual interaction, we use the seniority force in the T = 1 
channel, neglecting the neutron-proton correlations given 



the large value of T z = (N — Z)/2, and we calculate the 
strength with a given pairing window in the same approach 
as the one used by Moller and Nix [ll|. In the present 
study we choose to include in the pairing window all the 
single-particle states below the Fermi level and those ly- 
ing Ae = 6 MeV at most above the Fermi level. As in 
Ref. [l2j we have included a smooth cut-off as a function of 
the single-particle state energy with a diffuseness parameter 
/i = 0.2 MeV. 

We assume that GS nuclear shapes possess left-right re- 
flection and axial symmetries. Global studies [15 , Il4| indi- 
cate that close to the valley of /3-stability this is true except 
for a few small regions. Therefore the projection K of the 
total angular momentum of the nucleus on the z-axis, cho- 
sen to coincide with the intrinsic symmetry axis, and the 
parity 7r are good quantum numbers. In mean-field models 
this also holds for the single-particle states since the mean- 
field possesses these symmetries. For odd-mass nuclei, the 
time-reversal symmetry is broken but axial symmetry is 
assumed to be preserved. Depending on the intrinsic nu- 
clear deformation, different assumptions for the coupling 
schemes have to be considered. 

The nuclei for which the lowest-energy configuration is 
deformed are assumed to be rigid rotors. We can there- 
fore describe the coupling between the unpaired particle 
and the rotation of the even-even core in the rotor-plus- 
quasiparticle approximation using HFBCS single-particle 
states as in Rcf. [7j, and we assume that vibrations are 
"frozen" in the deformed well (zero phonon). To determine 
the lowest-energy quasiparticle in the HFBCS approach, we 
use the equal-filling approximation (see Ref. [151 ] and ref- 
erences quoted therein) to obtain a time-even state having 
the desired odd number of particles on average and then 
we create the lowest-energy quasiparticle on this state. 

For a spherical nucleus, the GS spin J and parity tt arc 
assumed to be those of the single nucleon, deduced from the 
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FIG. 1: (color online). Comparison with experimental data of the calculated GS spin and parity of odd-mass nuclei for which 
P2 ^ [0.01; 0.1] and V12 ^ 2 MeV throughout the nuclide chart within the HFBCS model using the SIII interaction. 



TABLE I: Agreement of GS spin and parity calculated with 
three different Skyrme interaction in the HFBCS approach and 
with the FRDM model with respect to Nubase 2003 data [2], 
for spherical ("Sph."), well-deformed nuclei ("Def.") and all of 
them ("Total"). For each Skyrme interaction, we indicate in 
parenthesis on the first line the percentage of agreement ob- 
tained when considering either of the two lowest quasiparticles 
and on the second line the number of nuclei for which the GS 
spin and parity of the second lowest quasiparticle agrees with 
the experimental data. 



Model Sph. Def. Total 

SIII 83.9% (90.8%) 40.5% (61.5%) 66.4% (79.0%) 
183(+15)/218 60(+31)/148 243(+46)/366 
SkM* 76.2% (89.2%) 37.5% (61.8%) 63.3% (80.0%) 
218(+37)/286 54(+35)/144 272(+72)/430 
SLy4 77.8% (85.8%) 39.3% (60.7%) 64.1% (77.6%) 
186(+19)/239 57(+32)/140 243(+51)/379 
FRDM 90.9% 43.1% 54.4% 
90/99 137/318 227/417 



quantum numbers (n, £, j, m with usual notation) of the 
last filled orbit, namely J = j and it = (— l) e . In our mean- 
field approach, this nucleoli occupies the first empty level 
of the self-consistent one-body potential above the lowest- 
energy levels occupied by the nucleons of the even-even 
core. 



In the FRDM model the single-particle states are ob- 
tained by diagonalization of the folded- Yukawa one-body 
Hamiltonian associated with the GS shape as explained in 
detail in Ref. (TB]. The ground state of the nucleus is ap- 
proximated by the Slater determinant built from the lowest- 
energy single-particle states (exhibiting the Kramer's de- 
generacy) . The spin and parity of an odd nucleus are there- 
fore those of the level occupied by the unpaired nuclcon, 
that is, the highest occupied level for the particle type in 
odd number [171 ] . This approximation is implemented for 
spherical as well as deformed nuclei. 

Within the HFBCS model, for computation time reasons, 
we search for the lowest-energy solution as a function of P2 
defined as in Ref. [IU and calculated exactly for the equiv- 
alent spheroid having the same quadrupole moment and 
mean square radius as the actual nucleus. The variational 
character of the HFBCS approach ensures that we obtain 
a minimum with respect to all the other shape degrees of 
freedom (compatible with the restrained symmetries). In 
principle we should search for the lowest local minimum in 
the whole multi-dimensional potential-energy surface, as in 
the FRDM approach. However this is a formidable task in 
a Hartree-Fock-like approach and not necessary in many 
cases. Indeed for several nuclei in various mass regions we 
have compared HFBCS higher multipole moments from the 
minimization in the $2 direction with FRDM ones obtained 
from a minimization in the full deformation space and found 
similar nuclear shapes. 

Since the potential-energy landscape is affected by the 
pairing correlations and since the strength of the seniority 
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FIG. 2: (color online). Same as Fig. [T]with the SkM* interaction. 



force to describe these correlations is not known a priori 
for a given nucleus, we have to proceed in an iterative 
way. Starting with the average set of pairing strengths 
for neutrons and protons, g„ = 17 MeV/(ll + N) and 
g p = 17 MeV/(ll + Z) respectively, we determine the local 
minima of the HFBCS deformation energy as a function 
of fa- Then we calculate the pairing strengths from the 
solution corresponding to the lowest minimum using the 
pairing model of Mollcr and Nix [ll[. With the new pair- 
ing strengths, we determine again the lowest local minimum 
and repeat this procedure until the P2 value of the GS min- 
imum converges to within 0.01. The pairing treatment in 
the FRDM approach is described in detail in Ref. [l3T ]. 

To solve the Hartree-Fock equations, we diagonalizc 
the one-body Hartree-Fock Hamiltonian in the cylindri- 
cal harmonic-oscillator basis. The truncation of the ba- 
sis and the approximate optimization of its parameters b 
and q are carried out as in Ref. Although the oscil- 

lator parameter b and the basis size parameter Nq scale 
as A" 1 / 6 and A 1 / 3 , respectively, we choose larger values 
of iVo as a function of A, namely values linearly inter- 
polated between iVo = 12 for the mass number A = 10 
and No = 18 for A = 260, and the optimal value of b for 
the spherical shape is taken to be the same for all nuclei, 
namely bo = 0.475fm _1 . We have checked that the GS 
radii, quadrupole moments, spins and parities of several 
nuclei across the nuclear chart are not sensitive to the basis 
size when we retain the above value of bo- 

We restrict the comparison of the HFBCS results with 
experimental data to nuclei for which the above models are 
expected to be valid. We therefore discard the nuclei for 



which 0.01 < P2 < 0.1 (interpreted to be soft) and those 
for which the energy difference V12 between the lowest two 
minima is less than 2 MeV (to avoid the possible ambigui- 
ties associated with shape coexistence). 

In Table|T]we show the corresponding percentage of agree- 
ment and the number of successful spin and parity calcu- 
lations with the SIII, SkM* and SLy4 interactions in the 
HFBCS model and with the FRDM model, for spherical 
and deformed nuclei. The same type of information is dis- 
played for each nucleus in the (N, Z) plane, separately for 
each model approach (three HFBCS and the FRDM calcu- 
lations) on Figs. [T] to [U 

In deformed nuclei, single-particle levels often lie very 
close together, which impairs our ability to correctly de- 
duce GS spins and parities. This is particularly true in 
heavy nuclei since the single-particle level density increases 
proportionally to A on average. Such an uncertainty in 
the determination of the GS spin and parity does not come 
only from the interactions used, but also from a possible re- 
ordering of nuclear levels in odd nuclei with respect to the 
order deduced from single-particle "bare" spectra caused 
by the coupling of single-particle degrees of freedom with 
other (collective) excitation modes. To estimate this uncer- 
tainty in the HFBCS calculations we consider not only the 
lowest-energy quasiparticle but the second lowest quasipar- 
ticle excitation as well when it lies at most 1 MeV above the 
lowest one. The gray (blue online) dots in Figs.[T]to[3]corre- 
spond to the nuclei for which the GS spin and parity of the 
second lowest quasiparticle state agree with the measured 
values, and the number of such cases is indicated in paren- 
thesis in Table [I] for each model and deformation category. 
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FIG. 3: (color online). Same as Fig. [T]with the SLy4 interaction. 
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FIG. 4: Same as Fig. [Q within the FRDM model. 
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In the first line for each model in Table [I] we also indicate in 
parenthesis the percentage of agreement when considering 
either of the two lowest quasiparticle excitations. 

As for the coupling with the core rotational mode in de- 
formed nuclei, the specific case of single-particle states with 
K = 1/2 has drawn our attention. Depending on the ac- 
tual value of the decoupling parameter a, one may have 
nuclear GS spins different from 1/2. This happens when 
a < — 1 or a > 4 and corresponds to about 30% of the cases 
where K = 1/2. In fact we have calculated a only for the 
lowest-energy quasiparticle when it has K = 1/2, so that 
the actual spin of the second lowest quasiparticle when it 
has K =1/2 has not been calculated. Therefore the num- 
ber of gray (blue online) dots in Figs. [1] to [3] is a lower limit. 
However this does not impair the conclusion drawn about 
the uncertainty arising from the high single-particle level 
density in deformed nuclei. 

As seen from Table [J and Figs. [1] to [3] the overall agree- 
ment is rather similar for all considered Skyrme forces. 
One may note, however, that the SIII parameterization is 
slightly better than the other two. In all four approaches 
the agreement for the spherical nuclei is much better than 
for deformed nuclei (partly because of the uncertainty men- 
tioned above for large single-particle level densities). For 
the spherical nuclei, the FRDM model is more successful 
than the HFBCS one (90% as compared to 80%, respec- 
tively) but with a much smaller set of nuclei, and yields 



similar results for deformed nuclei, taking only into account 
the lowest quasiparticle state. 

In the mass region A $J 100, including all spherical and 
deformed nuclei, the agreement is excellent in all HFBCS 
cases and slightly less good with FRDM (especially for the 
heaviest nuclei of this mass region). In deformed nuclei in 
or close to the rare-earth region, FRDM and SIII calcula- 
tions are equally good and yield slightly better results than 
the other two Skyrme forces, while in the actinidc region 
the best results are obtained from FRDM and SkM* calcu- 
lations (with a slight advantage for the former approach). 

To conclude this study, we should stress that given the 
global character of the study (of the order of 400 nuclei in- 
volved) , the assumptions made to connect calculated single- 
particle properties with observed spins and parities and 
the demanding character of such a comparison, the agree- 
ment can be deemed significant. This gives thus a rather 
good level of confidence in the predictive power of these ap- 
proaches when used in particular contexts. Improvements 
of this comparative study should include a more sophisti- 
cated treatment of the coupling of single-particle and col- 
lective degrees of freedom. 
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